Abstract 



Implementing the fanout gate by a Hamiltonian 

Stephen A. Fenner* 
University of South Carolina 

February 1, 2008 

co : 

o . 
o ■ 

(N ; 

. We show that, for even n, evolving n qubits according to a simple Hamiltonian can be used 

to exactly implement an (n + l)-qubit parity gate, which is equivalent in constant depth to 
an (n + l)-qubit fanout gate. We also observe that evolving the Hamiltonian for three qubits 
results in an inversion-on-three-way-equality gate, which together with single-qubit operations 
I . is universal for quantum computation. 
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■ 1 Introduction 

Let 7i be the Hilbert space of n + 1 qubits. The fanout operator F n+ \ : H — > H, depicted in 
^T) . Figure 1, copies the (classical) value of a single qubit to n other qubits. 

Fanout gates have been recently shown to be very powerful primitives for making shallow 
' quantum circuits [3, 4, 5, 6]. It has been shown that in the quantum realm, fanout, parity (see 

below), and Mod,j gates (for any q > 2) are all equivalent up to constant depth and polynomial 
size [3, 5]. That is, each gate above can be simulated exactly by a constant-depth, polynomial-size 
quantum circuit using any of the other gates above, together with standard one- and two-qubit gates 
(e.g., CNOT, H, and T). This is not true in the classical case. Furthermore, using fanout gates, 
in constant depth and polynomial size one can approximate n-qubit threshold gates, unbounded 
. \ AND (generalized Toffoli) gates and OR gates, sorting, arithmetical operations, phase estimation, 

^ ■ and the quantum Fourier transform [4, 6]. Since long quantum computations may be difficult to 

maintain due to decoherence, shallow quantum circuits may prove much more realistic, at least in 
the short term. 

On the negative side, fanout gates so far appear hard to implement. There is recent theoretical 
evidence that fanout gates cannot be simulated in small depth and small width, even if unbounded 
AND gates are allowed [2]. 1 All these results underscore how crucial fanout gates are to providing 
powerful, small-depth quantum computation. 

Rather than trying to implement fanout with a traditional small-depth quantum circuit, we 
can perhaps take another promising approach: evolve an n-qubit system according to a (hopefully 
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x Any circuit of depth d using AND and single-qubit gates to compute n-qubit fanout provably needs at least n/2 d 
ancillas. 
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Figure 1: Definition of the fanout gate. 



implementable) Hamiltonian. We show that a simple Hamiltonian, similar to one suggested recently 
by Chuang [1], does exactly this. 

Let YjZ be the n-qubit operator \ Ya=i ^i-> where Zi is the Z-gate acting on the ith qubit. If, 
for instance, each qubit is represented by a spin-^ particle, then is the observable representing 
the total spin angular momentum in the z-direction. Let J > be some constant in units of energy. 
We will show that for even n, the (n + l)-qubit fanout gate arises naturally by evolving the first n 
qubits through the Hamiltonian 

H n := J(EZ) 2 

twice — first for time |j then for time ^j- — together with a modest amount of additional processing. 
We also show that evolving the 3-qubit Hamiltonian J(T,Z) 2 results (modulo a global phase factor) 
in a 3-qubit "inversion on equality" gate /=, which maps \abc) to (— l) SabSac \abc) . Applying I = with 
one of the three qubits set to |1) results in a controlled Z-gate on the other two qubits, which is 
easily converted to a CNOT gate. 

Chuang's proposed Hamiltonian is closely related to H n . It is 

K n \= y j-ijZiZj, 

l<i<j<n 

where the Jjj are energy coefficients, and may be potentially realizable for certain combinations 
of the Jjj, depending on the physical arrangement of the qubits [1]. K n is the sum of pairwise 
interactions between the particles. In the special case where all the Jij are equal to J/2, we see 
that K n differs from H n by a multiple of the identity, and so evolving through H n and evolving 
through K n are equivalent up to an overall phase factor. 

We will show the 7= -gate in Section 2.2. In Section 2.3 we show the implementation of F n+ i 
for even n. In the sequel, we will assume for convenience that h = J/2 = 1. If X and Y are vectors 
or operators, we say, "X oc Y" to mean that X = e td Y for some real 9, that is, X = Y up to an 
overall phase factor. We use the same notation with individual components of X and Y, meaning 
that the phase factor is independent of which component we choose. If A is a set, we let || A|| denote 
the cardinality of A. 
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2 Main Results 



2.1 Evolving H n 

Since H n is represented in the computational basis by a diagonal matrix, it is particularly easy to see 
how it evolves in time. We only need to find the value of each diagonal element. Let 
be a vector of n bits. For 1 < i,j < n, we see that ZiZj\x) = (— 1) x *® x j \x), and thus ZiZj flips the 
sign iff Xi Xj. Suppose k of the xi are 1 and the rest (n — k) are 0. Then, 

ran 

(x\H n \x) = -££(-1)*^ 
i=i j=i 

= (l/2)(||{(i,j):x i =x J -}||-||{(t,j): a r i ^x i }||) 
= (l/2)(k 2 + {n-kf -2k{n-k)) 
= (1/2) (n 2 +Ak 2 -Ink) 
= n 2 /2 - 2k(n - k). 

The n 2 /2 term is independent of x, so up to addition of a multiple of /, the diagonal term is 
effectively —2k(n — k). We evolve H n for time t = 7r/4. Let U n := e~ iHnt where t = ir/4. We get 

(x\U n \x) oc e^ fc («- fc )/ 2 = (1) 
All the off-diagonal matrix elements of t/ n are zero. 

2.2 Implementing the J = -Gate 

For the n = 3 case, (1) yields 

f/ 3 ocdiag[l -1 -1 -1 -1 -1 -1 1 ] ocI=. 
No further gates are required. 

2.3 Implementing the i^+i-Gate 

Let n be even. is actually more closely related to the (n + l)-qubit (classical) parity gate, shown 
in Figure 2. We will show how to implement parity with U n . This suffices, because the equality 
given in Figure 3 was observed in [5, 3], where H is the Hadamard gate. 

Suppose for now that n = 2 (mod 4), and that x has k ones and n — k zeros. If k is even, then 
so is n — k, and thus k(n — k) = (mod 4). If A; is odd, then clearly, k(n — k) = 1 (mod 4). Thus 
we have, by (1), 

. _» rT ._, ( 1 if x has even parity, 

(x U n x) cc { . ., ^, AA ■/ 2 
[ i it x has odd parity. 

We start with the (n + l)-qubit basis state \x\ ■ ■ ■ x n ~irb), where x\, . . . , x n _i, r, b G {0, 1}. We 
prepare the nth qubit in the state 

2T|r> = (|0> + (-lHl»/>/2, 
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Figure 2: Definition of the parity gate. 
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Figure 3: Implementing fanout by a parity gate and Hadamard gates. 
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Figure 4: U n implementing the parity gate. 



|0> 



H 



Un 



Xi 



5+ 



H 



5+ 



^n— 1 



3?1 © • • • © %n—l 



Figure 5: A circuit similar to the parity gate. 



so that the current state of the first n qubits is \x\ 
state through U n , then by (2) the result is 



x n _i)(|0) + (-l) r \l))/V2. If we run this 



V2 



%1 ' ' ' %n—l 



)(^io)+^(-ini)), 



where p = (xi + ■ ■ ■ + x n -\) mod 2. The two states of the rath qubit corresponding to the two 
values of p are orthogonal, so we just need to rotate the nth qubit back to the computational basis. 
The nth qubit state is either +y or — y on the Bloch sphere, depending on p and r, so we can use 
HS^, where S is the phase gate diag [ 1 i ] , to rotate the y-axis to the z-axis. The final circuit 
is shown in Figure 4. When the CNOT gate is applied, its control qubit can be seen to be in 
the state i p \p(Br), unentangled with the other qubits. The rest of the circuit is then needed to 
uncompute the conditional phase factor V . We can implement C/ n ' by evolving — H n for time vr/4, 
or equivalently, by evolving H n for time 3ir/ 4, since = I. Of course, if we are willing to keep 
the phase factor, we can get by without this part of the circuit. In fact, we can get the parity-like 
gate of Figure 5 with only one use of U n . 

As we mentioned before, one gets a fanout gate by applying Hadamards on each qubit on both 
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Figure 6: U n implementing the fanout gate. 



sides of a parity gate. So starting with the circuit in Figure 4, after some simplification we get the 
circuit shown in Figure 6. If n = (mod 4), then we get a similar analysis of U n , except that (2) 
becomes 

. _. rT . f 1 if x has even parity, 
(x\U n \x) ot < ,, , J 

[ -i it i lias odd parity. 

It follows that we can swap U n with U n ^ in Figures 4, 5, and 6 above to maintain the equalities. 



3 Conclusions and Further Research 

A key point in our implementation is that the number of terms in the Hamiltonian H n is quadratic 
in n, which gives a quadratic term in the phase shift. We suspect there is also some way to get parity 
from the more general K n when the Jij are not all equal, provided there are still quadratically 
many terms. We also suspect that this will not work where there are fewer than quadratically many 
terms, for example, in the case where we just consider interactions between adjacent paricles in a 
ring, i.e., Jjj = J>0ifj = i + l (mod n), and Jij = otherwise. In this case, there are only 
linearly many terms in the Hamiltonian. 

A second Hamiltonian described by Chuang as potentially realizable in the lab [1] is 

L n = Ji,j{XiXj + YiYj + ZiZj), 

l<i<j<n 

which, when all the Jij = 2, differs from the squared total spin L 2 := (T,X) 2 + (T,Y) 2 + (T,Z) 2 by 
a multiple of /. We conjecture that L 2 can also be used to implement parity. 
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